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Abstract 



A phenomenological model for the inertial range scaling of passive-scalar 
turbulence is developed based on a bivariate log-Poisson model. An ana- 
lytical formula of the scaling exponent for three-dimensional passive-scalar 
turbulence is deduced. The predicted scaling exponents are compared with 

experimental measurements, showing good agreement. 
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The inertial range dynamics of a passive scalar advected by fluid turbulence is of great 
theoretical interest |l| . In particular, recent success @ in the deduction of the anomalous 
scaling exponent for the passive scalar field advected by a white and Gaussian velocity has 
inspired renewed enthusiasm for the problem ||. 

For the three dimensional turbulence diffusion, the advective velocity field is governed by 
the Navier-Stokes equation whose statistics is far from white and Gaussian. The extension of 
Kraichnan's theory [§] for this problem is intriguing and non-trivial. In fact it has even been 
a difficult task to experimentally measure the scaling exponents of the passive scalar. The 
existing experimental data do not yield convergent results • Direct numerical simulation 
(DNS) has been a useful alternative for studying scaling behaviors in fluid turbulence 
In a recent note, using the DNS data, we have tried to explain the discrepancies among 
experimental results 0. 

The existing phenomenological theories, such as the f3 model |9j and a bivariate log- 
Normal model fllQf , are based on fractal structures of turbulence ||, yielding analytical 
predictions of scaling exponents. Nevertheless, it has been noted that log-normal distribution 
leads to negative values of scaling exponents for high order passive scalar structure functions, 
contradicting existing numerical and experimental measurements. 

In this letter, we develop a phenomenological theory based on a bivariate log-Poisson 
model. The analytical formula for the scaling exponents is obtained. The scaling exponents 
from the theory are compared with those from experimental measurements, showing good 
agreement. 

According to Kolmogorov's refined similarity hypothesis (RSH), or K62 theory [jTT]] for 



the velocity field and similar RSH for the passive scalar field proposed by Obukhov |L2 

Au e = v u (te t ) l/ \ 
AT, = v^^N 1 / 2 , (1) 

where Aug = w(x + €) — u(x) and AT, = T(x + £) — T(x) are the longitudinal velocity 
increment and the scalar increment respectively, / is in the inertial range, v u and and vt are 
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random variables whose statistics only depend on the Reynolds number, eg and Ng are the 
locally averaged velocity dissipation and scalar dissipation, respectively [|J. Equation (|]) 
has been verified by experiments and is in general believed to be correct [[HJ. Suppose that 
the p — th order velocity and passive scalar structure functions satisfy the scaling relations 
in the inertial range: S p (£) = ((Aug) p ) ~ fa, S^{£) = ((ATg)' p ) ~ £ z ", then to deduce scaling 
exponents from (jl|), a probability density function (PDF) of velocity field and a joint PDF 
for velocity and passive scalar must be modeled. Here (•) denotes an ensemble and £ p and 
z p are the p — th order scaling exponents for the velocity and scalar structure functions, 
respectively. The log-normal assumption in K62 theory [[□]] has been challenged recently by 
She and Leveque JE1] who argued that a log-Poisson PDF for e r leads to a better description 



of the statistics of e r . A hierarchy model based on the log-Poisson PDF was derived for 
the inertial range scaling exponents, showing good agreement with numerical simulations |7j 



and experiments [15 



The distribution function of a random variable x satisfying Poisson distribution can 



be written as ||16|| : p(x — i) — Pi — e a a*/i!, where a is the variance and the mean for 
the variable. The generating function of pi has the form: P(s) = Y^iPiS 1 = e~ a+as . The 
generating function for a bivariate Poisson distribution p it j = p(x — i,y = j) can be defined 

as 



a-L—a2—b+aisi+a2S2+bs\S2 



where ax, a 2 are variances for x and y, respectively, and b is a constant representing the 
correlation between x and y. Let us assume that the velocity dissipation and passive scalar 
dissipation have the hierarchy structure relation: eg 1 = Wg x g 2 tg^ and Ng x = V^N^, where C\ 
and £ 2 are two length scales; Wg Y g 2 and Vg^ are multiplicative factors depending on £ 1 and 
i 2 . For Ei and l 2 within the inertial range, the multiplicative factors can be written in the 
forms of Wg x g 2 = (£i/£ 2 ) h /3 x and Vg x g 2 = (£i / £ 2 ) hT (3% , where (3, (3t, h and h T are constants 
to be determined later; (X, Y) are stochastic variables. In the following, we assume that 
(X, Y) follow the bivariate Poisson distribution pij. 
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If there is a scaling, then for £ in the inertial range, (efiVj) ~ £ T ( p ' g \ where t va is 
the scaling exponent of power order p and q. Using the bivariate Poisson distribution 
assumption, we have: 

It is easy to recognize that the summation in the right hand side of the above equation 
is nothing but the generating function for the bivariate Poisson distribution. It is then 
straightforward to write a formula for r(p, q): 

rip, q) = -hp -h T q + (a + a T + b- a/3 p - a T (3 q T - bp p /ln{£ x / '£ 2 ), (2) 

where constants a, cit and b are functions of £± and £ 2 . The obvious physical constrains 
r(l,0) = r(0, 1) = result the following relations: a = [h/(l — /3)}ln(£ 1 /£ 2 — b) and ax = 
[/&t/(1 — l3T)}ln{£\ / £ 2 — b). If we further assume the correlation between ei and Ng satisfying 
power law scaling in the inertial range, i.e. b = 7b/n(£i/£ 2 ) and 75 is a constant, we obtain 
the scaling exponent relation: 

r(p, q) =r p + r q + r b (p, q), (3) 

where 

T p =-h P +[h/{l-j3)](l-p), 

and 

T q = -h T q+[h T /(l-P T )](l-ft), 
are scaling exponents for moments (ef) and (iV|) respectively, and 

r 6 (p,g) = 7b (l-/3 p -^ + W) 
is the contribution from the cross correlation. There have been several discussions about 



how to determine the model constants in t va |14],|17| . The parameter (5 depends on the 
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constant h and the intermittency parameter // = T2 H17| , i.e. (5 = 1 — fi/h. It has been 
argued that the co-dimension, C, for the most singular structure of the velocity dissipation 



TM is related to (3 and ft.: 



C = h/(l-f3) = h 2 /fi. 



(4) 




-0.3 - 



FIG. 1. Theoretical scaling exponents (p < 0) for the velocity structure function field as a 
function of p (x), compared with numerical measurements (□) 



The original hierarchy structure model by She-Leveque assumes that the high intensity dissipa- 
tion region tends to form tube-like structures which implies: C = 2. Using fi = 2/9 |J9], we obtain 
/i = ±2/3. The sign associated with h is due to the square operation in (|3|). If we further assume 
z p grows slower than exponential as p tends to ±00, the positive (or negative) h has to be chosen 
for the positive (or negative) power p in @). The original SL model with h = 2/3 corresponds to 
the range of positive powers, yielding (3 = h = 2/3. As shown in (|l|) for Sp with p > 0, only the 
negative power p for r p is involved. Hence, for p < 0, we choose h = —2/3, yielding j3 = 4/3 and: 



Tp = -p + 2[l-(4/3f]. 



(5) 



Using RSH relation in ([!]), we have 



C p = -p + 2[l-(4/3y/ 3 ]. 



(6) 



In a recent paper [18|, we have compared DNS data with the theoretical prediction using h = 2/3 
for (p when —1 < p < 0. In Fig. 1, we show again this comparison with (||) for h = —2/3. 
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The similar procedure can be carried out for the passive scalar dissipation. For q > 0, only the 
positive q in r q has contribution. Consequently, a positive Ht has to be chosen. Using the similar 
arguments as for the velocity field, we finally obtain Ct = hx/0- — Pr) = h T /nT, where ht is the 
intermittency parameter for the passive scalar dissipation Ng. It is noted that the high intensity 
region for the passive scalar dissipation N tends to form sheet-like structures J8|, implying the 
co-dimension Ct = 1. We then obtain Ht = yfj^r and j3r = 1 — y / /^r- From numerical || and 
experimental || data, [it ~ 1/4, we have: 

T q = - l -q+[l-{l/2)% (7) 

From (|5|) and (|7|), r(p,q) can be written as follows: 

r{p, q) = 3 + \p - 2(~f - ig - (1/2)" + 76 [1 - (^f - (1/2)" + (|y»(l/2)«]. (8) 

Using (||) and (^), we have the following analytical formula for z p : 

z P = 3P + t(— I , |) 

3 - ±p - 2(3/4)f/ 6 - (1/2)p/ 2 
+ 76 [1 - (3/4)f/ 6 - (1/2)p/ 2 + (3/4)f/ 6 (l/2)f/ 2 ]. (9) 

In the above equation, 7^ is not determined. At this moment, there is no direct measurement of 
7b and we treat it as a free parameter. In Fig. 2a we plot z p as functions of p for 7 ^ = 0,0.1 
and 0.2. The scaling exponents from the experimental measurement carried out by Antonia et. 
al. H have also been included. A good agreement between theory and measurement for 7^ ~ 
is seen, indicating that the correlation between q and N% is weak. In fact, the variation of the 
scaling exponents is not drastic for 7& varying from to 0.2. In Fig. 2b, we demonstrate the 
model dependence on co-dimension Ct- We notice that when Ct is between 0.75 and 1, the model 
results agree with the experiment. As Ct increases further, the model prediction deviates from 
the experimental data. It has been found that the co-dimension of passive scalar field is around 
0.6 jilj. It would be useful to directly measure Ct, the co-dimension of passive scalar dissipation 
using experimental or numerical simulation data. 
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The evidence of scaling exponents for 75 = matches well with those from experimental 

measurements might indicate real physics. As a matter of fact, from (|]) we notice that for positive 

and large p values, the major contribution of e e p ^ 6 to ST comes from small amplitude events whose 

structures are blob-like, whereas the major contribution of N^ 2 to Sj comes from large amplitude 

events whose structures are sheet-like. The former is often associated with weak stretching regions 

while the latter is associated with strong stretching regions. We suspect that the weak spatial 

correlation of these regions might be the direct reason why the nearly-independent log-Poisson 

with 7^ ~ gives a good prediction of the scaling exponents. We should mention that (H|) fails to 

be valid if power p is too large. In fact, when p is larger than 36 and 75 equals to zero, z p starts be 

negative. Up to now, there is very little understanding about the statistics of very high moments, 

due to the difficulties in both experiments or numerical simulations. Consequently, there is no solid 

justification for using the RSH and Hierarchy model in such high power range. 
2 
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FIG. 2. Theoretical scaling exponents z p for the passive scalar field as functions of p compared 
with experimental measurements 1|, showing the model dependence on: (a) 7^ while the Ct equals 
1; (b) Ct while 76 equals 0. 
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